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Abstract 

Background: Improvement in muscle strength is an important aim for the rehabilitation 
of hemiplegia patients. Presently, the rehabilitation prescription depends on the 
evaluation results of muscle strength, which are routinely estimated by experienced 
physicians and therefore not finely quantitative. Widely-used quantification methods for 
disability, such as Barthel Index (Bl) and motor component of Functional Independent 
Measure (M-FIM), yet have limitations in their application, since both of them 
differentiated disability better in lower than higher disability, and they are subjective and 
recorded in wide scales. In this paper, to explore finely quantitative measures for 
evaluation of muscle strength level (MSL), we start with the study on quantified 
electromyography (EMG) and sonomyography (SMG) features of tibialis anterior (TA) 
muscles among hemiplegia patients. 

Methods: 12 hemiplegia subjects volunteered to perform several sets of plantar-flexion 
movements in the study, and their EMG signals and SMG signals were recorded on TA 
independently to avoid interference. EMG data were filtered and then the root-mean- 
square (RMS) was computed. SMG signals, specifically speaking, the muscle thickness of 
TA, were manually measured by two experienced operators using ultrasonography. 
Reproducibility of the SMG assessment on TA between operators was evaluated by 
non-parametric test (independent sample T test). Possible relationship between muscle 
thickness changes (TC) of TA and muscle strength level of hemiplegia patients was 
estimated. 

Results: Mean of EMG RMS between subjects is found linearly correlated with MSL 
(R^ = 0.903). And mean of TA muscle TC amplitudes is also linearly correlated with MSL 
among dysfunctional legs (R^ = 0.949). Moreover, rectified TC amplitudes (dysfunctional 
leg/ healthy leg, DLHL) and rectified EMG signals (DLHL) are found in linear correlation 
with MSL, with R^ = 0.756 and R^ = 0.676 respectively. Meanwhile, the preliminary results 
demonstrate that patients' peak values of TC are generally proportional to their personal 
EMG peak values in 12 dysfunctional legs and 12 healthy legs (R^ = 0.521). 

Conclusions: It's concluded that SMG could be a promising option to quantitatively 
estimate MSL for hemiplegia patients during rehabilitation besides EMG. However, after 
this exploratory study, they should be further investigated on a larger number of subjects. 

Keywords: Sonomyography, Tibialis anterior muscle. Surface electromyography, Muscle 
strength level 
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Background 

Hemiparesis, or muscular weakness of the upper and lower legs contralateral to the 
brain lesion, is the most frequent sign following a stroke or intracerebral haemorrhage 
(ICH) [1], The severity of paresis post stroke is an important consideration since it is a 
common manifestation of stroke that relates to functional capabilities, discharge destin- 
ation and mortality [2]. Therefore, improvement in muscle function is an important 
aim during the rehabilitation of hemiplegia patients. The rehabilitation prescription de- 
pends on the evaluation results, which are routinely estimated by experienced 
physicians. 

At present, there are many clinical evaluation methods of disability measure for 
stroke patients with hemiplegia, such as Fugle-Meyer method [3], Bobath level method. 
National Institute of Health (NIH) stroke scale [4], Scandinavian Scale [5], Brunnstrom 
stage method and so on [6]. The measures of motor performance that they incorporate 
do not challenge muscles sufficiently to indicate their strength accurately [7]. Re- 
searchers or practitioners may have difficulties in interpreting the clinical meaning of 
the scores or changes in scores. Moreover, different physicians may make different 
evaluation decisions because the evaluation process is not finely quantitative. Addition- 
ally, the physicians can't make individualized prescription due to the lack of detailed 
movement information [7]. 

To quantify disability in patients, activities of daily living (ADL) (Barthel Index, BI, 
and motor component of Functional Independent Measure, M-FIM) and categorical 
disability measures (Modified Rankin Scale, MRS) are also used [8]. Although used glo- 
bally, the BI and M-FIM scales have limitations in their applications [8]. Yu et al. [7] 
presents Brunnstrom stage automatic evaluation for stroke patients using extreme 
learning machine. The real time movement data is collected by using 3-axis acceler- 
ometer sensors which are placed in the geometric center of arms. 

Muscle strength testing has been used to monitor status and recovery of hemiplegia 
patients [9]. Considerable evidences suggest that the ability to perform a physical task 
is determined by a threshold level of muscular strength and endurance [10,11]. A de- 
cline in functional status is determined at least in part by muscle strength, flexibility, 
range of motion, physical fitness, and body composition [10,12]. Although older co- 
horts have the highest risk of developing disability, the association between muscular 
strength and endurance and the subsequent prevalence of functional limitations found 
in [13] indicates that this relationship persists even among middle-aged adults. The 
findings suggest that maintenance of strength throughout the lifespan may reduce the 
prevalence of functional limitations. 

A wddely-used tool to quantify disability in hemiplegia patients is EMG. Many studies 
have been previously reported on the relationship between the surface electromyogram 
(EMG) and muscle force [14-16], length [17], wrist joint torque [18], etc. Most studies relat- 
ing the surface EMG to muscular force suggest that the amplitude of the EMG should in- 
crease proportionately with the square root of the tension, rather than linearly, if motor 
units fire independently of one another [19]. Some reports found that mean force-EMG re- 
lationships were highly non-linear but similar in shape for different cats and different speeds 
of locomotion [16]. But some authors have concluded that the magnitude of the EMG sig- 
nal is directly proportional to muscle strength for isometric and/or isotonic contractions 
with constant speed for various muscles. Direct experiments have shown a linear 
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relationship [20,21]. The force and integrated EMG relationship averaged across subjects 
showed almost the same linear correlation in spite of different electrode locations [22]. No 
linear correlation was observed between median power frequency of EMG and the knee ex- 
tension force [22]. Although surface EMG is used to quantify muscle activity, the relation- 
ship between force and surface EMG during voluntary contractions is not fully understood 
[23]. EMG has been shown to be a useful method of evaluating patients with stroke and 
other neurological disorders [24,25], but factors that prevent the direct quantification of 
muscle force from EMG signal include cross-talk [26], variations in the location of the re- 
cording electrodes [27] and the involvement of synergistic muscles in force generation. And 
its difficult for surface EMG to detect the deep muscles non-invasively, due to the fact that 
the deep muscle EMG may be attenuated more and mixed by the superficial muscle EMG 
when reaching the skin surface [28]. 

Ultrasound is capable of visualizing normal and pathological muscle movements [29-42]. 
Ultrasound appeared to be even more sensitive in detecting fasciculations compared to 
EMG and clinical observations, because it can visualize a large muscle area and deeper lo- 
cated muscles [29,41]. Neuromuscular disorders give rise to structural muscle changes that 
can be visualized with ultrasound: atrophy can be objectified by measuring muscle thick- 
ness, while infiltration of fat and fibrous tissue increase muscle echo intensity [29]. There- 
fore, besides EMG, there is a newly prevailing means to quantify disability in patients, 
named as sonomyography (SMG). SMG [30], defined as the muscle thickness change de- 
tected from the real-time ultrasound images was related to the morphological change of the 
muscle during its contraction. SMG has been employed for the control of prosthesis [30], 
the assessment of muscle fatigue [31], isometric muscle contraction [32], dynamic muscle 
contraction [33], one dimensional SMG tracking the guided patterns of wrist extension 
[28], the relationship between torque and rectus femoris muscle morphological change in 
healthy adults [34], estimation and visualization of longitudinal muscle motion [35] and so 
on. However, there is a lack of research using SMG to investigate the relationship between 
thickness change (TC) of tibialis anterior (TA) muscles and muscle strength level (MSL) 
hemiplegia patients. 

The purpose of this study was to investigate the relationship between TA muscle s 
EMG /SMG features and MSL during plantar-flexion contraction. The relationship be- 
tween EMG and SMG features was also studied. 

Methods 

Subjects 

Some patients who had cerebral surgery and came to the physical therapy rehabilitation 
center within one month volunteered to participate in the experiment in Guangzhou 
Zhujiang Hospital, Guangdong, China. The study was approved by the institutional review 
board of Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences. 

Enrolment was fulfilled the following inclusion criteria: (1) stable vital signs within 1 
month; (2) cerebral infarction (CI) and intracerebral haemorrhage were verified by com- 
puted tomography angiography (CTA), magnetic resonance angiography (MRA) or mag- 
netic resonance imaging (MRI); (3) patient was at the stage for post CI/ICH rehabilitation 
admission in the hospital. The exclusion criteria included: (1) cognitive disorder which 
made patient incapably follow the guidance of examiner; (2) severe motor dysfunction 
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required contact guard; (3) paraplegia & quadriplegia; (4) course of CI/ICH over 3 weeks; 
(5) motor dysfunction diagnosis before the onset of CI/ICH; (6) patient required ankle-foot 
orthotic; (7) any physical complaints impeded the test, such as muscle pain. 

Consent forms were obtained from the patients and his/her authorized legal repre- 
sentative. November-December 2012, a group of 27 patients were in the physical ther- 
apy rehabilitation center for post CI/ICH rehabilitation. 12 patients fulfilled conditions 
and participated in the study (9 male, 3 female, age 52.2 ± 11.1 years). The number of 
dysfunctional left legs of patients is 6, and the number of dysfunctional right leg is 6. 
The general information of patients in different muscle strength level (MSL) is listed in 
Table 1. 

Experimental procedure 

Patient seated on an anchored wheel chair performed ankle joint dorsiflexion following 
the direction of examiner. The movement was completed under the command of 
examiner. Patient performed dorsiflexion with no force restriction. The excessive knee 
and hip joint movement were prohibited during plantar-flexion. Practice was carried 
out in ahead of test for patient knowing the start position, lift and lower action in 
dorsiflexion. The qualified movement was ensured in practice. Any undesired move- 
ments including overactivated/depressed activity and ankle eversion/inversion in dorsi- 
flexion was corrected by examiner. If the undesired movement was caused by CI/ICH, 
it would be cancelled without any correction. 

a) Start position: Patient kept the sole of foot on the floor and hold leg perpendicular 
with the floor (Figure 1). 

b) Lift action: First activity (begin with the tone of start') to lift the sole of foot away 
from the ground while heel was inactivated on the floor. During this action, 
patients perform the maximal voluntary contraction (MVC) of the feet. 

c) Lower action: Second activity (begin with the tone stop') of the foot back to the 
ground. 

At least three consecutive trials should be completed. A two minutes rest was given 
after practice to avoid muscle fatigue in the test. Root mean square (RMS) of EMG was 
used to check the muscle activity. The qualified signal was recorded by the examiner. 

EMG assessment 

Subject was seated in a chair with the knee flexed at 90 degrees and the ankle at neutral 
position. The surface EMG signals were recorded from the TA of right leg (channel 1), 

Table 1 Patients medical record in Guangzhou Zhujiang Hospital 



Muscle Strength Level (MSL) Number of subjects Number of dysfunctional legs (left/right) 



1 


1 


1/0 


2 


3 


2/1 


3 


2 


1/1 


4 


5 


1/4 


5 


13 


1/0 
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Figure 1 Patient's EMG signal acquirement and ultrasound acquisition. 

J 



the TA of left leg (channel 2), the gastrocnemius of right leg (channel 3) and gastrocne- 
mius of left leg (channel 4) by four pairs of surface electrodes (Figure 1). Each electrode 
(BioNomadix 2-Channel EMG system, BIOPAC Systems, Inc., Goleta, CA, USA) was 
separated from the other by 20 mm. Sampling rate was set at 1000 samples per second 
(1 KHz). The fixed gain of 2000 and 10 Hz high pass, 500 Hz low pass and notch filter 
off were preset for filter option. Noise voltage was 1.5 [iV rms with bandwidth of 
1.0 Hz to 500 Hz. The EMG data were exported as a MATLAB file and analyzed for 
peak value. EMG signal was quantified from RMS with interval of 0.5 s. 

Maneuvers were carried out three times each, with the patient sat in a comfortable 
examination desk. Subjects were asked to initiate voluntary motor task movements 
under commands from the examiner for the duration of the flexion movement. 



SMG assessment 

SMG assessment was not synchronized with EMG to avoid interference. Linear array 
ultrasound (US) transducer (Sonostar Technologies Co., Guangzhou, Guangdong, 
China) was placed perpendicular to the TA between the site of EMG and elec- 
trodes for transverse imaging measurement (Figure 1). The transducer was held prop- 
erly by examiner to ensure consistent images in the test. Activity was following the 
command so that examiner could hold the position of probe on target. Besides, ad- 
equate coupling gel was used on transducer to guarantee sound wave well travel 
through the skin. Frequency of US transducer was at 7.5 MHz with 70 mm depth de- 
tection. Anatomical cross-sectional area (ACSA) was shown in Figure 2. Focus, contrast 
and brightness were adjusted for different imaging. In the trial, 128 frames within 10 
seconds were collected for SMG analysis. 

Images were manually measured in matlab (The Math Works Inc., Natick, MA, USA). 
Image was cropped to a suitable size and analyzed in the sequence of frames. The TA 
ACSA was delineated by using solid white line and separated by muscle fascia boundary 
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and the other tissues' architecture in the images. Two examiners had been trained and 
experienced on muscle architecture finding in US image independently and were 
blinded from the profiles of patients. The difference between maximum thickness and 
minimum thickness of TA muscle thickness AB as showed in Figure 2 were measured 
twice. AB was recorded with pixel numbers of 6 per (mm). Thickness change (TC) be- 
tween TA muscle in activity period (AP) and TA muscle in rest period (RB) was calcu- 
lated by (1) of 

TC = (ABniaximum~ABniinimum)/6 (l) 

AB ranging was admitted to visually differentiate the AP and RB by manual plot of 
continuous US images. The AP and RB periods were selected by examiner through ob- 
servation that there were no paired points with the increasing and decreasing tendency 
indicating another AB ranging period in AP or RB. 
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Statistical analysis 

Statistical analysis was carried out in the SPSS version 19.0 software. Non-parametric 
test was applied to analyze SMG ACSA. Reproducibility of the SMG assessment on TC 
between two operators was evaluated by non-parametric test (independent sample T 
test). Two experienced examiners assessed each image two times in two different days. 

Results 

EMG analysis 

Excel 2007 was used for the EMG analysis and SMG analysis. The EMG signal of pa- 
tient was recorded during three-times plantar- flexion movements (Figure 3). It was 
found that mean of patients' EMG amplitudes in each scale of MSL was proportional 
to MSL (Figure 4). The relationship for EMG RMS versus MSL was best- fitted with lin- 
ear model (R^ = 0.903). 

SMG analysis 

To investigate the significant difference between two examiners' TC results, two group 
samples from different examiners were assessed by independent sample T test result. 
Statistical significance was accepted at the 5% level. There was no significant difference 
between two operators' results (p > 0.1). 

To study the relationship between mean of TA TC in different hemiparetic legs and 
MSL, we computed the mean values of TA TC in each scale of MSL. It was observed in 
Figure 5 that mean of TC amplitude of patients was in good correlation with the MSL 
(R^ = 0.949). Mean of TC amplitude of hemiplegia legs increased almost linearly as the 
MSL increased. 

At the same time, variances of EMG and SMG for all subjects in the same MSL were 
found large. Figure 6 and Figure 7 show boxplots of EMG RMS with peak value in the 
same MSL and the mean of TC amplitudes in the same MSL respectively. 

In order to explore the relationship between mean of EMG RMS and mean of TC 
amplitude, a chart of the mean of EMG amplitudes and corresponding TC amplitudes 
is shown in Figure 8. A very high linearly correlation could be concluded between these 
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Figure 4 Mean of different patients' EMG amplitudes versus muscle strength level. 



two quantified parameters (R = 0.957). Furthermore, to study individual patient's rela- 
tionship between EMG signal and SMG signal, all EMG amplitudes and TC amplitudes 
(R^ = 0.521) in 12 dysfunctional legs and 12 healthy legs were shown in Figure 9. There 
is a nearly linear relationship between the EMG amphtude and TC amplitude for each 
patient. EMG amplitudes and TC amplitudes in 12 dysfunctional legs have a higher lin- 
ear relationship (R^ = 0.686). Notice that it indicated that each patient' inherent correl- 
ation between the EMG signal and TC signal in dysfunctional leg. 

Considering the individual difference among hemiplegia patients, a method was to be 
explored using the value of TC result in dysfunctional leg normalized by TC in healthy 
leg of same subject. The intercept of the linear regression of percentage (dysfunctional 
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leg/healthy leg, DLHL) was used to calculate the representative value. Rectified TC am- 
plitudes (DLHL) were proportional to MSL as shown in Figure 10. The relationship of 
rectified TC amplitudes (DLHL) and MSL was fitted with a linear model (R^ = 0.756). 
Moreover, mean of rectified EMG RMS amplitudes (DLHL) versus MSL was studied. 
There was a similar relationship between rectified EMG RMS amplitudes (DLHL) and 
MSL (R^ = 0.676) in Pi gure 11. Figure 12a and Figure 12b show boxplot of rectified TC 
amplitudes (DLHL) and rectified EMG amplitude (DLHL), respectively. 

Based on the results from Figures 4, 5, 6, 7, 8 we presented reasonable hypothesis 1 
that TC of patient in dysfunctional leg was smaller than that of patient in healthy leg. 
Given the hypothesis 1 exists, decisions are made to differentiate the dysfunctional leg 
and healthy leg with the TC value. Prediction results of the evaluation theme are shown 
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in Table 2, The accuracy is 83.3% and 91.7% for operator 1 and operator 2 s result, 
respectively. 



Discussion 

The report on TC during plantar-flexion among hemiplegia subjects is limited, there- 
fore in this study, we engaged two experienced operators to process the US images 
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Figure 9 EMG amplitudes and corresponding TC amplitudes in 12 dysfunctional legs and 12 
healthy legs. 
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manually. Though there were no statistical differences between their outputs, and we 
have taken the averaged values from both operators as TC, it should be noted the 
inter-operators difference is not totally trivial, 0.52 ± 1.39 mm. 

Generally speaking, as shown in Figures 4 and 5, results of both EMG, as a popular 
means of muscle study, and SMG (TC of TA), are in strong correlation to the clinically 
prevailing MSL, which provided proofs to the assumption that the MSL can also be 
reflected muscle bio-electric activities and morphological changes. 

However, it should be noted that in this experiments, large inter-subjects variances 
were observed for both EMG and SMG. For EMG, it has previously reported that the 
variance of EMG may provide additional information that is "lost" when using only an 
average EMG amplitude as EMG variance may indicate how the average muscle activity 
is occurring [43]. And fortunately, EMG and SMG themselves demonstrated good cor- 
relation, as shown in Figure 8, though they represented two distinct genres of signals 
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Figure 12 Boxplot of percentage (DLHL) (a) of patients' TC amplitudes, (b) of patients' 
EMG amplitudes. 



generated during the muscle contractions. Furthermore, for individual patients, EMG 
amplitudes and corresponding TC amplitudes also had a linear relationship (as shown 
in Figure 9) in totally 12 dysfunctional and 12 healthy legs. 

If we excluded the possible measurement errors, a major cause of the large inter- 
subjects variances therefore can be the variances of execution of the exercises. Specific- 
ally speaking, the differences of subjects in terms of both cognitive status and level of 
willing to wholly complete the exercises, as hardly controllable in an objective way, 
could lead to the large inter-subjects variances of EMG or SMG. In other words, if a 
subject does not try his/her best of strength, a variance will be contributed to the 'true' 
relationship between this scale of MSL and EMG or SMG, but less possible the rela- 
tionship between EMG and SMG themselves. 

Another concern is the 'true' inter-subjects variance for each MSL. Hence as stated 
in the method section, we also evaluated the DLHL. A highly significant relationship 
between rectified TC amplitude (DLHL) and MSL was observed (Figure 10). And a 
highly significant relationship between rectified EMG amplitude (DLHL) and MSL was 
observed (Figure 11). 

The results of the experiments to differentiate dysfunctional leg and healthy leg using 
the TC showed a very good performance, as shown in Table 2, thus provided prelimin- 
ary supports to hypothesis 1. 

Previous research had shown that deficits in EMG area were not observed in all par- 
etic muscles. The type of contraction influenced the degree of deficits demonstrated in 
the paretic muscles, with greater preservation of force during eccentric contractions 
[44]. The EMG of the paretic TA could be used as a control signal for the beginning of 
electrical stimulation, through the detection of changes in their RMS values during the 
hemiparetic gait [45] and TA EMG data could be used in sleep/wake estimation [46]. 
The reason we choose TA as the object muscle in this study is, TA is independent 



Table 2 Prediction results of the evaluation theme based on hypothesis 1 



Operator 


The number of correct_predicted 
subjects 


The number of wrong_predicted 
subjects 


Accuracy 

(%) 


#1 


10 


2 


83.3 


#2 


11 


1 


91.7 
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relatively and less possibly affected by other muscles during the plantar-flexion move- 
ment. Yet large variance existed in this study for the TA EMG signals as shown before, 
and the reasons remain unclear and require further investigations. 

Conclusions 

In an exploratory study among the hemiplegia patients, both mean of EMG RMS in 
each MSL scale for TA and mean of TA TC in each MSL scale during plantar-flexion 
were strongly correlated with MSL. There was a good correlation between mean of 
EMG amplitudes and TC amplitudes as well. There was a nearly linear relationship be- 
tween the EMG amplitude and TC amplitude for each patient among dysfunctional 
legs. Results from the current study provided preliminary insights to the quantitative 
relationship between EMG/SMG and MSL as well as EMG and SMG themselves, in an 
example muscle of TA producing the plantar-flexion exercises among hemiplegia pa- 
tients. Thus it suggests that SMG could provide an alternative perspective to under- 
stand the muscle contraction activity morphologically among hemiplegia patients, 
along with the more well-reported EMG. Further investigations on how such morpho- 
logical information could be applied to evaluation of muscle function, yet require ex- 
periments on a larger number of subjects with further clustering according to factors 
such as age, gender or pathological origins. 
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